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1. INTRODUCTION 

An electron liquid on a uniform neutralizing background (EL) is used as the 
reference system in most realistic calculations of electronic structure in condensed- 
matter physics [1]. At zero temperature there are only two relevant parameters 



for a disorder-free EL in the absence of quantizing magnetic fields and spin-orbital 
coupling: (i) the usual Wigner-Seitz density parameter = (ymaDft^)"^^^, cis being 
the Bohr radius in vacuum (in a medium — > = fi^«;/(mbe^) with R and mb 
appropriate dielectric constant and bare band mass respectively); and (ii) the degree 
of spin polarization ( = \n^ — n^l/naD- Here Uu is the average density of particles 
with spin a =t, J, and 77,20 = + ni is the total average density. 

Understanding the many-body aspects of this model has attracted continued 
interest for many decades [2-6]. The EL, unlike systems of classical particles, behaves 
like an ideal paramagnetic gas at high density (rg <^ 1) and like a solid at low 
density [7] (r^ ^ 1). In the intermediate density regime, which is relevant in three 
dimensions (3D) to conduction electrons in simple metals and in two dimensions 
(2D) to electrons in an inversion layer of a Si metal-oxide-semiconductor field-effect 
transistor (MOSFET) or in an AlGaAs/GaAs quantum well, perturbative techniques 
are clearly not effective owing to the lack of a small expansion parameter. One has to 
take recourse to approximate semi-analytical methods, a number of which have been 
reviewed in Refs. [3,4], or to Quantum Monte Carlo (QMC) simulation methods [8- 
18], 

Among the methods designed to deal with the intermediate density regime, of 
particular interest for its physical appeal and elegance is Landau's phenomenological 
theory [19]. The basic idea of Landau's theory is that the low-lying excitations of a 
system of interacting Fermions with repulsive interactions can be constructed starting 
from the low-lying states of a noninteractig Fermi gas by adiabatically switching-on 
the interaction between particles. This procedure allows to establish a one-to-one 
correspondence between the eigenstates of the ideal system and the approximate 
eigenstates of the interacting one. Landau called such single-particle excitations of 
an interacting Fermi-liquid "quasiparticles" (QP's). He wrote the excitation energy 
of the Fermi-liquid -E'[np] as a functional of the QP distribution function Up in terms 
of the isolated quasiparticle energy £p and of the QP-QP interaction function /p,p'. 
The latter can in turn be used to obtain various physical properties of the system, 
such as the compressibility and the spin-susceptibility. 

One of the implications of Landau's theory of Fermi-liquids is the fact that the 
QP mass m* is renormalized by electron-electron interactions [2] : m* 7^ m, m being 
the bare electron mass. In a translationally-invariant system the current jp carried by 
a single excited QP of momentum p is controlled only by the bare mass, jp = p/m. 
On the other hand, the QP group velocity Vp is instead defined by Vp = Vp£^p. In 
an isotropic system Vp is parallel to p and the relation Vp = p/m* defines the QP 
effective mass. Thus jp ^ Vp, the reason being that due to interactions the moving 
QP tends to drag part of the electronic medium along with it producing an extra 
current [2,3]. 

The QP effective mass is a measurable quantity. The most direct way to deter- 
mine m* would be a measurement of the low-temperature heat capacity Cy(T). It 
is in fact remarkable [2,3] that electron-electron interaction effects enter Cv{T) only 
through m*: the Landau interaction functions /p,p' are not explicitly invoked [2,3]. 
These type of experiments are exceedingly challenging and have not yet been realized 
with high precision (sec Ref. [38] below and references therein). An alternative tool 
to access experimentally the QP effective mass (and other Fermi-liquid parameters) 
is to analyze quantum Shubnikov-de Haas oscillations of the magnetoresistance. Mo- 



tivated by a large number of recent magnetotransport studies [20-25] of this type 
we have revisited [26] the problem of the microscopic calculation of m* in a para- 
magnetic {i.e. C = 0) 2-D EL. Our systematic study is based on a generalized GW 
approximation which makes use of the many-body local fields and takes advantage 
of the results of the most recent QMC calculations of the static charge- and spin- 
response of the 2D EL [27]. We report extensive calculations for the many-body 
effective mass enhancement m* /m over a broad range of Tg values. In this respect we 
critically examine the relative merits of the on-shell approximation, commonly used 
in weak-coupling situations, versus the actual self-consistent solution of the Dyson 
equation. We show that already for ~ 3 and higher, a solution of the Dyson equa- 
tion proves here necessary in order to obtain a well behaved effective mass. Finally 
we also show that our theoretical results for a quasi-2D EL, free of any adjustable 
fitting parameters, are in good qualitative agreement with some recent measurements 
in a GaAs/AlGaAs heterostructure. 



2. MANY-BODY EFFECTIVE MASS ENHANCEMENT 

The aim of this Section is to present a brief summary of the theory we have used 
for the retarded QP self-energy Ei.et(k, cu) of a paramagnetic 2D EL (that we have 
summarized in Eqs. (3) and (7) below) from which we have calculated m*/m. The 
formal justification of Eqs. (3) and (7), which essentially rests on both diagrammatic 
perturbation theory and on the so-called renormalized Hamiltonian approach [28], 
can be found in the original works [29,30], in Ref. [3], and in Ref. [26]. 

To fix the notation we start by introducing the Hamiltonian for a 2D EL confined 
to an area S, 

7i = XI ^k4,aCk,a + ^ I] 4i+q,ai cL-q.a^Ck^.a^ Cki ,ai • (1) 

k,cT q^O ki,CTik2,cT2 

Here Cj[ ^ and Ck,^ are fermionic creation and annihilation operators which satisfy 
canonical anticommutation relations, £k = /i^k^/(2m) is the single-particle energy, 
and Vq = 2'Ke^ jq is the 2D Fourier transform of the bare Coulomb interaction 
e^/r. For later purposes we introduce the Fermi wave number kp — (27rn2D)^'^^ = 
^/2/{rsaB)i the Fermi energy sp = fT?k\/{2m) and the quantity ^k = £k — £f- 
The retarded QP self-energy Sret (k, a;) is written as the sum of two terms, 

Sret (k, Lo) = Esx (k, oj) + Ecu (k, a;) (2) 

where the first term is called "screened-exchange" (SX) and the second term is called 
"Coulomb-hole" (CH). The frequency cj is measured from sp/h. 

The SX contribution is given by 

Ssx(k, u)^- [ e(-a+q//^) . (3) 

J (27r)2 £(q,a;-Ck+q/^) 

Here Q{x) is the step function and £(q, a;) is a screening function originating from 
effective Kukkonen-Overhauser interactions [31], 

= 1 + [1 - 6"+ (q)] ' Xc (q, a;) + 3 (q) Xs (q, a;) • (4) 



In Eq. (4) the charge-charge and spin-spin response functions xc(q)'^) and Xsi^j^) 
are determined by the spin-symmetric and spin-antisymmetric static [32] local-field 
factors G+(q) and G-(q) [27], 



1 -'f;q[l -G'+(q)]xo(q,^) 

and 

. X xo(q,^) 

Xs(q,ct') = — „ . . — 7 r, (6) 

l + 'f;qG'_(q) Xo(q,a;) 

Xo(q, being the Lindhard-Stern response function of a noninteracting 2D elec- 
tron gas [34]. Note that Esx(k, a;) is just an ordinary exchange- like self-energy built 
from the Kukkoncn-Overhauser effective interactions instead of bare Coulomb inter- 
actions, which would lead to the frequency-independent Hartree-Fock self-energy first 
calculated for the 2D EL by Chaplik [35]. 

The CH contribution to the retarded self-energy is given by 



ScH k, = - / — •Wq / jr 77— TT, (7) 

where 5 is a positive infinitesimal. The real part of the retarded self-energy is readily 
obtained from Eqs. (3) and (7) with the result 

5ReSret(k,a;) = - j ^ ^q5Re[£-i(q, a; - a+qA)] ©(-^k+qA) 

Once the real part of the QP self-energy is known, the QP excitation energy 5£^Qp(k), 
which is the QP energy measured from the chemical potential /i of the interacting 
EL, can be calculated by solving self-consistently the Dyson equation 

where 3?eE^^(k, a;) = 3?eEret(k, a;) — SretC^F, 0). For later purposes we introduce at 
this point the so-called on-shell approximation (OSA). This amounts to approximat- 
ing the QP excitation energy by calculating 3f?eE^^(k, u) in Eq. (9) at the frequency 
u) — ^\<i/h corresponding to the single-particle energy, that is 

5£:Qp(k) Ck + 5ReS^,,(k,a;)L^^^/, . (10) 

The effective mass m* can be calculated from the QP excitation energy by means 
of the relationship 

m* Ti^kF dk k=kF 



As remarked above, 5SQp{k) may be calculated either by solving self-consistently 
the Dyson equation (9) or by using the OSA in Eq. (10). In what follows the identity 



will be used, uj{k) being an arbitrary function of k. 

Using Eqs. (11) and (12) with u;{k) — SSQp{k)/h we find that the effective mass 
calculated within the Dyson scheme is given by 



D 



m 1 + {m/h^kF) dkUei:f^,{k,uj%^ 



(13) 



The renormalization constant Z that measures the discontinuity of the momentum 
distribution at A; = A;f is given by 

^= i-h-^d^^ei:UKu^%=,^^^^^- ^^^^ 

The normal Fermi-liquid assumption, < Z < 1, implies dJR,eT,^^^{k,uj)\^^^^ < 
0. Thus we see that the effective mass rn^ can diverge at a finite value of 
by one of two mechanisms: (i) the partial derivative of with respect to uj^ 
9t^3?eEj^^(/c, a;) |^_^^ going to minus infinity at some finite value of [36]; (ii) 

the partial derivative of Sj^^ with respect to /c, 9fc3fJeS^^(A;, a;)|^_^^ going to 

—h^kF/m at some finite value of r^. 

Neither possibility is realized in our calculation: the first is barred a priori by 
the fact that the analytic expression for the frequency derivative of is always 
finite at finite rs', the second is found a posteriori not to occur since the momentum 
derivative of is positive up to the largest rs considered (see below). 

On the other hand, using Eqs. (11) and (12) with u}{k) = ^k/^ we find that the 
effective mass TOosa within the OSA is given by 



OSA 



1 + {m/h^kF) d,^eEl,{k,oj)\ + (1 - Z-i) 



(15) 



Of course, Eq. (15) is a valid approximation to the effective mass in the weak coupling 
limit, as can be seen by expanding Eq. (13) for small values of E^^: however its 
application becomes problematic at large values of rs- In particular, we see that 
because Z decreases monotonically with increasing r^, there must necessarily be a 
critical value of for which the denominator of Eq. (15) vanishes and rriQgj^ diverges. 
A recent paper by Zhang and Das Sarma [37] infers from this fact a true divergence 
of the effective mass within the RPA. In our view, however, this must be considered 
an artifact of Eq. (15). The unphysical character of the behavior mQSA ^ oo is 
revealed by the fact that the divergence is driven by a negative but finite value of 
9cj3?eE^^(/ci;', 0), whereas we know, from the general analysis given above, that a 
genuine divergence would have to be driven either by an infinite di^^eTj^^^{kF, 0) or 



by a negative 9fc3?eE^^(fci?, 0) becoming equal to —h^kp/'m. We conclude that there 
is no evidence, within the present theory, for a physically relevant divergence of the 
effective mass. 



3. NUMERICAL RESULTS 

We turn now to a presentation of our main numerical results for m*/m. In all 
figures the labels "RPA", and refer to three possible choices for 

the local-field factors: "RPA" refers to the case in which local-field factors are not 
included, to the case in which the antisymmetric spin-spin local field is set to 

zero {i.e. spin-density fluctuations are not allowed), and finally ^^G^&lG^" refers to 
the full theory including both charge- and spin-density fluctuations. 

In Fig. 1 we show our numerical results for and TOqsa- The effective mass 
enhancement is substantially smaller in the Dyson-equation calculation than in the 
OSA, the reason being that a large cancellation occurs between numerator and de- 
nominator in Eq. (13). In both calculations the combined effect of charge and spin 
fluctuations is to enhance the effective mass over the RPA result, whereas the oppo- 
site effect is found if only charge fluctuations are included - a manifestly incorrect 
result that neglects the spinorial nature of the electron. For completeness we have 
also included in Fig. 1 the variational QMC results of Kwon et al. [11]. The reader 
should bear in mind that the effective mass is not a ground-state property and thus 
its evaluation by the QMC technique is quite delicate, as it involves the construc- 
tion of excited states. There clearly is quantitative disagreement between our "best" 
theoretical results (the "G_|_&G_/D" predictions) and the QMC data. 

In Fig. 2 we show the behavior of the two terms in the denominator of Eq. (15) 
as functions of r^. This flgure clearly shows how a spurious divergence can arise 
in mosA- instance, within the RPA the denominator in Eq. (15) has a zero at 
Tg — 15.5 (see the inset in Fig. 2). Our numerical evidence, within the three theories 
we have studied, is that indeed (i) dJRellf^^{kF,^) is negative as it should for a 
normal Fermi-liquid, and monotonically increasing in absolute value as a function of 
r^; and (ii) 9fc3?eSj|^(/cF, 0) is positive and monotonically increasing too. Within the 
theory outlined in Ref. [26] , which uses as a key ingredient the Kukkonen-Overhauser 
effective screening function in Eq. (4) , the effect of a charge-only local fleld is to shift 
this divergence to higher values of r^, while the opposite occurs upon including both 
charge and spin fluctuations. For instance, within the "G_|_&G_/OSA" theory the 
divergence occurs near = 5. 

We now pass to illustrate how our theory compares with experiments. A full 
analysis of the published data for the effective mass of carriers in Si-MOSFET's [20,21] 
would require a more complete theoretical study, mainly to account for the two- valley 
nature of the material. We will focus here instead on the experimental results of Tan 
et al. [38] in a GaAs/AlGaAs hetero junction-insulated gate field-effect transistor 
(HIGFET) of exceedingly high quality. A quantitative comparison between theory 
and experiment would also require a refined treatment of a series of effects such as 
those due to (i) the detailed band-structure of the host semiconductor, (ii) disorder, 
and (iii) finite temperature [39]. Note also that the Shubnikov-de Haas measure- 
ments of Tan et al. [38] are performed in a small but obviously finite magnetic field 




2 4 6 8 10 

Figure 1. Effective mass enhancement as a function of Vg for < < 10. The 
inset shows an enlargement of the results for < 1. The lines show the results from 
Eq. (13), while the symbols (except for the dots) are from Eq. (15). The QMC data 
(dots) are from Ref. [11]. 



B and, in general, the ground-state of an EL in a small finite B field [40] can be 
profoundly different from that at 5 = 0. In our zero-field calculations we have: (i) 
included band-structure effects only through the GaAs band mass mb = 0.067 m; 
(ii) neglected possible effects due to disorder because the concentration of back- 
ground impurities rii in the HIGFET used in Ref. [38] has been estimated [41] to 
be ~ 5 X 10^^ cm~'^ which is a very low number; and (iii) neglected thermal ef- 
fects even though in Ref. [38] the temperature of the dilution refrigerator was kept 
relatively high (100 ^ T < 400 mK) to avoid the quantum Hall regime in which the 
Shubnikov-de Haas oscillations become non-sinusoidal. A simple inspection shows 
however that, assuming that for the highest temperature T = 400 mK the EL is in 
thermal equilibrium with the refrigerator, the ratio of thermal to Fermi energy is 
quite small even at the lowest densities, e.g. kBT/sp ~ 0.03 for Vg = 6. 

We thus restrict our analysis solely to the effect of finite sample thickness, by dis- 
cussing how a softened Coulomb potential modifies m* against the strictly-2D results 
shown in Fig. 1. The expectation is that the QP effective mass will be noticeably 
smaller when a softened Coulomb interaction is at work. 

We have thus recalculated m* after renormalizing the bare Coulomb potential by 
means of a form factor to take into account the finite width of the EL in the HIGFET 
used in Ref. [38] . The appropriate renormalized potential is given by = VqJ^{qd) / k, 
where 
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Figure 2. Illustrating the divergence of the effective mass within the OSA. The three 
curves starting from unity at = refer to the quantity l+m^~^fc^^9fc5?eE^^ {kp: 0), 
and the other three curves to h~^\duj'SieTif-^^{kF,0)\ = — 1. The intersection of 
two lines with the same line-style in the two sets of curves corresponds to a zero in 
the denominator of Eq. (15) and thus to a divergence in mQSA- "^^^ inset shows this 
divergence occurring within the RPA at ~ 15.5. 

with d = [^^Ksc/(487rmbe^n*)]^/^ representing an effective width of the quasi-2il> 
EL [42]. Here Kins = 10.9 and Ksc = 12.9 are the dielectric constants of the insulator 
and of the space charge layer, h is their average and n* = ndepi + 11?T'2d/32, the 
depletion layer charge density naepi being essentially zero (see Ref. [18] of Ref. [41]) 
in the experiments of Ref. [38] . Note that the renormalized potential does not con- 
tain any adjustable fitting parameter. The results that we obtain with the softened 
potential are shown in Fig. 3 against the experimental results of Tan et al. [38] . 
A caveat to keep in mind is that we have used the same local-field factors as for a 
zero-thickness 2D EL [27] in the lack of a better choice. Thus the results labeled by 
"G-l-" and in Fig. 3 contain the effect of finite thickness only through the 

renormalization of the Coulomb potential. We believe that the explicit dependence 
of the local fields on the finite width of the 2D EL should not change the results of 
Fig. 3 in a substantial manner. 

As it appears from Fig. 3, a quite satisfactory qualitative agreement exists be- 
tween theory and experiment, considering the oversimplified model we have used for 
the quasi-2il> EL in the heterojunction of Ref. [38]. Nevertheless, in detail the dis- 
crepancies are quite substantial. In particular both the RPA and "G'-(-&G'_" theories, 
which give rather similar results, underestimate the considerable effective mass en- 
hancement at strong-coupling, i.e. for > 5. A comparison in the weak-coupling 
regime r^. < 1 would be very helpful but unfortunately is not possible due to the 
lackness of experimental data. 



1.5 
1.4 
1.3 

% 1.1k 
1 

0.9 
0.8 



RPA/D 
G+&G_/D 



RPA/OSA □ 
G+&G_/OSA ♦ 
Tan et al. i — •- 



I 



n LJ AT 








Figure 3. Effective mass enliancement for a 2D EL confined in a GaAs/AlGaAs 
triangular quantum well of the type used in Refs. [24] and [38]. The theoretical results 
are shown in the same notation of Fig. 1. The experimental results (Tan et al.) are 
from Ref. [38]. 



4. CONCLUSIONS 

In summary, we have revisited the problem of the microscopic calculation of the 
many-body effective mass enhancement in a 2D EL. We have performed a systematic 
study based on the many-body local-fields theory, taking advantage of the results of 
the most recent QMC calculations of the static charge and spin response of the EL 
expressed through static local-field factors. We have presented results for the effective 
mass enhancement over a wide range of electron densities. In this respect we have 
critically examined the merits of the OSA versus the Dyson-equation calculation. 
Depending on the local-field factors, the OSA predicts a spurious divergence of the 
effective mass at strong coupling and a solution of the Dyson equation is therefore 
necessary in order to obtain the correct value of the effective mass within Fermi-liquid 
theory. The comparison with the experimental data of Ref. [38] shows good qualita- 
tive agreement but substantial discrepancies especially at strong coupling calling for 
further theoretical work and computer simulations. 
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